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ABSTRACT 



Arguments for including in the VO large grids of syn- 
thetic spectra fully covering the HR diagram are pre- 
sented. One obvious need is that of population synthesis 
at all redshifts. Theoretical spectra have also the power to 
predict peculiar behaviors, that could remain unnoticed in 
observed spectra, or lead to erroneous conclusions. One 
interesting example is given of the Ca II H & K lines in 
extremely metal-poor stars. In carbon-rich atmospheres, 
the H & K lines become much weaker relative to the 
continuum, which will lead to an underestimate of the 
metallicity. What actually happens is a displacement of 
the continuum, obvious in absolute-flux synthetic spectra, 
but not visible in observed continuum-normalized spec- 
tra. Libraries of synthetic spectra are described, as well as 
the codes used to compute them. A few remarks are made 
on issues with making these libraries available through 
the VO, as well as on the necessary input data, i.e. line 
lists and model atmospheres. 

Key words: Star: atmospheres; Stars: chemically pecu- 
har; Technique: spectroscopic; Virtual Observatory. 



WHY SYNTHETIC SPECTROSCOPY IN THE 
VO? 



spectra are necessary to extract stellar parameters from 
observations (Tcs, [Fe/H], etc). 

Various kinds of spectra can be computed, corresponding 
to the usage that is intended: population synthesis may re- 
quire low resolution, full spectral coverage, or conversely 
a medium resolution, shorter spectral coverage (e.g. H 
band). High resolution, over a short spectral interval is 
required to prepai-e the GAIA RVS (Ca II IR triplet re- 
gion). Exploitation of VLT-UVES spectra requires very 
high resolution over the whole optical and near-IR do- 
main. In all cases, the variation of many parameters is a 
key to the superiority of synthetic spectra over observa- 
tions (Toff, gravity, [Fe/H], [a/Fe], C, N, O abundances, 
etc). This may represent a large number of spectra if a 
significant portion of the HR diagram is to be covered (at 
least lO'' to 10* spectra). The computation time of such 
an amount of spectra is not a big issue if the model at- 
mospheres used are classical (in the sense 1-D, LTE, hy- 
drostatic). If one at least of these simplifying hypotheses 
is lifted the computing time may become prohibitively 
large. For some regions of the HR diagram, or if a great 
precision is wanted, this proves necessary (e.g NLTE and 
wind hydrodynamics for OB stars, non-LTE and 3-D hy- 
drodynamical models of convection for better than 5% 
accuracy in abundances in solar type stars). What is al- 
ways necessary is a good set of physical input data (line 
lists, partition functions, etc). 



Synthetic spectra have several advantages over observed 
spectra: they naturally provide absolute fluxes (and in- 
tensities!), as well as normalization to the continuum, 
and they can be computed for whatever stellar parameters 
(even rare or unobserved types!). In principle, a database 
of synthetic spectra can cover the whole HR diagram 
at all metallicities, and non-standard chemical composi- 
tions, which is invaluable for stellar population synthesis. 
Synthetic spectra may be used for the preparation of in- 
struments, e.g. GAIA, or to understand characteristics 
of sub-samples in large surveys (using e.g. color-color 
diagrams or low-resolution spectra). They are necessary 
to quickly classify spectra in large surveys (e.g. micro- 
lensing searches) with non-standard wavelength ranges 
at various spectral resolutions. And of course synthetic 



2. WHAT CAN BE DONE WITH SYNTHETIC 
SPECTRA: FOUR ILLUSTRATIONS 



2.1. impact of stellar parameters on spectra 

I chose here to illustrate this on some unusual extremely 
metal-poor carbon-rich stars. Carbon-enhanced stars 
make up around 20% of very metal-poor stars and have 
become a boiling field of research in recent years (see the 
review of [Beers & Christlieb | ( |2005l l). The model atmo- 
spheres used for the calculation of the spectra shown on 
Fig. [T] |2l an d [3] we re computed with the MARCS code 
j^stafsson et al. | [T975) |Plez et al. 1992 ; Gustafsson] 
[etal. ,,2003) . The Fi gures speak for themselves. Similar 
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Figure 1. Effect of T^g variations on the optical and 
NIR spectrum of extremely metal-poor giant model stars 
(MARCS models). Note the disappearance of Hydro- 
gen lines (and jumps) and the strenghtening of molecular 
lines ( e.g. CH around 4300A with decreasing tempera- 
ture. All spectra on an absolute flux scale at the stellar 
surface. 



4000 6000 8000 lo' 1.2x\0^ 

wavelength (A) 

Figure 3. Same as Fig. [7]/or carbon abundance varia- 
tions. Molecular lines become progressively prominent at 
larger [C/Fe] ratios: CH around 4300A, or C2 at larger 
wavelengths. Note also the distortion of the UV contin- 
uum. 

observed spectra are used to search for metal-poor stars 
with specific characteristics of temperature, metallicity or 
carbon enhancement in low-resolution spectroscopic sur- 
veys Uke the SDSS-SEGUE (Sloan Extens ion for Galac- 
tic Understanding and Exploration, see e.g. |Sivarani et al.| 
[] ( 2005[ )). Also, synthetic photometry can be directly 
computed from these spectra and used for similar pur- 
poses. One particular interest is the calibration of color- 
temperature relations for extremely metal-poor stars, es- 
pecially of non-standard composition, that are only ob- 
served in small numbers, if at all. 




[Fe/H] = -3 [alpha/Fe]=04 
[Fe/H]=-5 [alpha/Fe]=04 5000K logg = 2.0 
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Figure 2. Same as Fig. |7J/or variations of [Fe/H]. Note 
the weakening of the Ca II triplet around 8500A. 



2.2. Effect of carbon enhancement on the Ca II H & 
K lines 

Fig. |4] shows the near-UV model spectra of extremely 
metal-poor main sequence stars with T'off=5000K. The 
Ca II lines appear weaker with higher C/Fe ratio, which 
gives the impression at first glance that the Ca abundance 
is lower As many searches for metal-poor stars have re- 
lied upon H & K strength (in particular the many can 
didates extracted from the HK survey of Beers et al. 
( 1985 1), it is important to understand this H & K weaken 



ing to be able to disentangle Ca and C abundance effects. 
What happens is understood by looking at the opacity 
in the continuum and line forming region in the mod- 
els. In the [C/Fe]=0 model atmosphere, the dominant 
continuum opacity source around 3500A is H^, whereas 
in the model with [C/Fe]=2, the CH continuous absorp- 
tion increases the opacity by a factor 10 to 100. This 
strong difference in the continuous opacity deeply affects 
the continuum flux, whereas line cores, with much larger 
opacities, remain mostly unchanged. This can be seen in 
Fig. [5] (absolute flux showing the change in continuum 
level), as well as in Fig. [3] (distortion of the continuum 
in the near-UV). The effect is only seen because the flux 
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Figure 4. Effect of carbon abundance variations on the 
Ca II H & K lines in an extremely metal-poor giant model 
star (MARCS models). The Ca abundance is the same in 
all models. Flux normalized to the continuum, as would 
be the case for observed spectra. 
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Figure 5. Same as fig. |4] but on an absolute flux scale. 
Note that the continuum level is very different in the four 
models, but the Ca II line cores are unaffected. The Ca 
abundance is the same, only the equivalent width of the 
lines is much smaller in the C-rich spectra. 



is plotted on an absolute scale, allowing comparison be- 
tween spectra of stars with different abundances. When 
observing at the telescope, the situation would look as in 
Fig.|4] and therefore look as a Ca underabundance, not a 
C overabundance! Only synthetic spectra can reveal such 
an effect. 



2.3. Other illustrations 



Further illustration is given by the possibility to predict 
theoretical equivalent widths of, e.g., Balmer lines as a 
function of Tcff in hot stars, useful in population synthe- 

(12005 ). A last exam- 



sis, as shown by Martins et al. 



pie is given by the wor k of [Massey et al. "| (|2005 ) and 



Levesque et al. (2005 2006 1, using synthetic MARCS 



spectra and observed spectro-photometry of red super- 
giant stars to derive simultaneously the Tcff of the stars 
and the reddening on the line of sight. The reddening law 
affecting the slope of the spectrum, does not change the 
strength of the TiO bands relative to the local continuum. 
As this strength is very sensitive to temperature, the Tcff 
can be derived independently of the reddening, that fol- 
lows from an adjustment of the slope of the theoretical 
spectrum on the observed spectro-photometry. This work 
has lead to a revision of the temperature scale of red su- 
pergiants. Another result is that the reddening towards 
red supergiants is larger than in the surrounding OB as- 
sociations, which shows that there is a circumstellar dust 
contribution to the reddening. The slope of the redden- 
ing indicates that the dust grains are of larger size than 
their interstellar counterparts. Such a work would not be 
possible without synthetic spectra, as virtually all red su- 
pergiants are reddened by unknown amounts, presumably 
by grains of various sizes. 



3. SOME CHALLENGES FOR 
SPECTROSCOPY IN THE VO 



SYNTHETIC 



Some of the challenges ahead are, .e.g., (i) the analysis 
of huge amounts of spectra from GAIA, SDSS -SEGUE, 
RAVE, and multifibre spectrographs, esp. from ELTs, 
each of which represents millions to billions of spectra, 
or (ii) to synthesize stellar populations at high redshifts, 
low metallicity, in various wavelength regions. For that, 
we need very large grids of synthetic spectra, and fast, 
reliable and accurate algorithms to do it automatically. 
Two approaches may be chosen: to use large libraries of 
spectra from the VO to compare to observations with VO 
tools or to compute specific spectra for each observation 
, either locally with resident software and computer, but 
using data and other tools from the VO, or with codes run- 
ning on a server inside the VO. We will further comment 
on these two philosophies below. 
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Figure 6. Sampled SED from a MARCS model compared 
to a high resolution spectrum computed with the same 
model and almost the same line data. The flux at each 
sampled wavelength is tlie same as that computed at high 
resolution, but it gives no clue to how the flux varies be- 
tween points. A local integration of an equivalent width 
will be erroneous, but on large scale the average flux level 
is correct. 



4. LIBRARIES OF SYNTHETIC SPECTRA 



A number of libraries already exist on the web, although 
not always integrated in the VO. Spectral Energy Distri- 
butions (SED) are available for the various brands of 1-D, 
s tatic, L TE models ATLAS ( jKurucz |[T979l|Castem eTaT] 
H [T997]l, MARCS (marcs.astro.uu.se), and PHOENIX 
dHauschildt et~ar~l [19%! [Hauschil dt et~al~~| [1999}. 
Note that in some cases the SEDs are actually sampled 
spectra, as the models are usually computed with the 
opacity sampling algorithm, a Monte-Carlo approach to 
the radiative field calculation. All radiative quantities 
are computed at a large number of wavelengths, but 
nothing is known in between. The difference between a 
high resolution spectrum, and a spectrum sampled at a 
resolution of R ^ A/ A A = 20000 is shown in Fig|6] It 
is clear that the sampled flux cannot be compared to an 
observed spectrum at the same R = 20000 resolution. 
This can be done only at a degraded resolution. A sam- 
pled spectrum will globally provide the right fluxes at a 
much lower resolution, but not locally at high resolution. 
Detailed high-resolution spectra are only available for 
subsets of the existing grids. There are too many to 
cite them all, but most can be found (together with ob- 
served libraries) on the excellent site of D. M. Gutierrez 
(www.ucm.es/info/Astrof/invest/actividad/spectra.html). 
A problem still to be solved for many of these libraries 
is to properly include them in the VO (see one approach 
in the contribution by |Rodrigo Blanco ( 2007) in this 
volume). 



most of the HR diagram, at all metallicities. It will also 
provide SEDs at lower resolution. Regular updates are 
planned. 



4.1. Issues to keep in mind 

Computing such large grids of models and spectra may 
become time-consuming. CPU time is of course needed, 
but also scientist time to assemble the input data, set up 
the calculations, and check that the output makes sense. 
Making the data available on the web, other than as a 
tar-file on an anonymous ftp, and especially if it has to 
be VO compliant may also become a burden. Here, help 
from our friends of the VO is critical. The benefits of 
getting ones models easily accessible on the web, are 
twofold: (i) the models get used and one hopefully gets 
acknowledged (the concern is that one has to find a way 
that when data is made available through the VO, and is 
passed from one tool to another to finally become part 
of some end product, proper referencing to the original 
source is ensured), (ii) as importantly, the use of the mod- 
els towards different purposes, sometimes not anticipated 
by their producers, leads to a thorough tests of these mod- 
els and their limits. However, users must remember that 
model spectra may be flawed, because of inappropriate 
physics (e.g. LTE O star atmospheres, ID static models 
for convective 3D atmospheres), because of inaccurate 
input data (e.g. line positions and strengths), or because 
of missing data (e.g. cool stars without TiO). It is impor- 
tant to understand the limits of a given brand of models 
and to refrain from using them to far off limits. Differ- 
ences will always appear at some level between observa- 
tions and calculations, but the problems are often not easy 
to ascribe to either models or observations. Observers 
tend to blame the models, and vice-versa... One way to 
estimate systematic errors in models may be to compare 
models computed with different codes and input data sets. 
While not always possible this is certainly very instruc- 
tive. Illustrations can be found in [Martins et al. | ( |2005| l. 
Differences tend to increase for more extreme spectral 
types, and progress needs to be made particularly on car- 
bon stars, and very cool stars, although successes like the 
analysis of very cool AGB star spectra with Teg down to 
2700 K by Garcia-Hernandez et al. ( 2007| l demonstrate 
the level of sophistication of modem model atmospheres 
and spectra. 



One library of synthetic spectra that is just 
about to be made available is the POLLUX database (see 
the contribution by |Palacios et al. ] (|2007) in this vol- 
ume). It will be directly VO compliant and includes high 
resolution synthetic spectra (300-1200 nm), covering 



5. SYNTHETIC SPECTROSCOPY AT HOME OR 
SPECTRA ON DEMAND FROM VO ? 



To address the challenges outlined above, one approach 
would be to implement in the VO interfaces to synthetic 
spectroscopy codes that could be run on demand by the 
user for specific combinations of stellar parameters. This 
would be very flexible, but prohibitively demanding on 
resources due to the huge amount of computations that 
could be required (new calculation for each slightly dif- 
ferent stellar parameter or abundance). A better way of 
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doing is to provide in a simple and standard way large 
grids of pre-computed spectra, at very high resolution (it 
is easy to degrade them to lower resolutions), with only 
the main stellar parameters varied (Toff, gravity, Fe/H, a- 
elements/Fe, and possibly some CNO variations). These 
grids can be interpolated to produce population synthe- 
ses. All present needs in population synthesis can pre- 
sumably be fulfilled with these basic grids. Similarly, 
to analyse large numbers of observed spectra, the calcu- 
lated grids can be interpolated to feed automatic analysis 
codes, resulting in first estimates of the main stellar pa- 
rameters. Refinements can be made locally using a spec- 
trum synthesis code, model atmospheres, and spectral 
line lists. The automatic tools for stellar parameters deter- 
mination are certainly something one would want to see 
integrated into the VO. The development of these codes 
has followed various lines: minimum distance methods, 
that are limited to few parameters, neural networks that 
may be difficult to teach, and also new solutions that seem 
fast a nd efficient Uke MATISSE ( |Recio-Blanco et al. 1 
|2006t . 



5.1. Synthetic spectroscopy codes (LTE) 

If the user wants to compute spectra with specific stel- 
lar parameters, we have seen that the only viable solu- 
tion is a local implementation of a synthetic spectroscopy 
code. However, ingredients to compute synthetic spec- 
tra may be found on the web, in the VO: stellar model 
atmospheres, and atomic and molecular line lists (see 
below). A number of synthetic spectroscopy codes are 
available from their authors. I will here only present 
the 1-D LTE codes that are most widely used (multi- 
D codes, or polarized radiation codes are usually not 
freely available; some non-LTE codes can be found on 
the web). SYNTHE was developed by R. Kurucz and 
is adapted to the calculation of spectra for the ATLAS 
models. It is available on line either on Kurucz page (ku- 
rucz. harvard.edu), or better, with updates and documen- 
tation, on wwwuseroat.ts.astro.it/atmos. MOOG was ini- 
tially written by C. Sneden, who continues to update it. 
The present version is pure-LTE (no continuum scattering 
included), and intended for FGK stars. It is available with 
documentation on verdi.as.utexas.edu/moog.html. SME 
(spectroscopy made easy) is developed by N. Piskunov, 
and is further described in his contribution in this vol- 
ume. It allows in particular abundance stratifications. 
It is available upon request. TURBOSPECTRUM was 
developed by B. Plez from an older code of the Upp- 
sala group (B. Gustafsson). The assumptions are similar 
to the MARCS code and the input physics is the same. 
It is optimized for the computation of cool star spec- 
tra (FGKMSC) and can compute large chunks of spectra 
with millions of spectral lines. It is available upon re- 
quest, without documentation. Some issues have to be 
kept in mind when using off-the-shelf codes. Making a 
code available to a wide audience is not a simple task, 
if one wants to avoid its use in conditions it was not in- 
tended for Many warnings and conditional tests have to 
be implemented (this is usually not done, and the result 



will be that the code may run and produce an erroneous 
output without any warning). Writing a documentation 
is a tedious task. It may be incomplete, and is often not 
read thoroughly. Ideally, there should be continued con- 
tact between users and providers of codes in order to help 
improve, correct, debug and document them. 



5.2. Spectral line databases 

A main ingredient for the computation of syn- 
thetic spectra is the line lists. There are many 
databases on the web. There are also lists 
of databases, always incomplete, but useful: 
www.cfa.harvard.edu/amdata/ampdata/amdata.shtml. 
The Kurucz home page is an enormous resource, but not 
always well documented with regards to, e.g., what input 
data was used to construct the lists. An excellent atomic 
line database is the VALD (jPiskunov et al. , 1995 ; Kupka] 



etal. 1999 Piskunov 2007 1. A third release is planned 
for the Fall 2007, including for the first time molecular 
lines (TiO). Molecular line lists are scattered on various 
sites: NASA- Ames, University College London, ... A 
good source of information is the stellar atmosphere code 
descriptions, as their conceptors strive to use the best and 
most complete data. The contribution by E. Roueff in 
this volume provides also a few more references. 



6. SUMMARY AND A FEW QUESTIONS 



I have shown that the setup of large libraries of synthetic 
spectra is needed in the VO and quite straightforward to 
implement if only a few stellar parameters are allowed to 
vary. Of course the data model needs to be finalized, but 
the discussions we had at the conference will certainly 
speed up the process. For example, calculated spectra 
cannot be queried with a position on sky, which is the 
usual way to query observations... On the other hand 
one would like to be able to query synthetic spectra with 
an abundance pattern, or a temperature. Such libraries 
will allow to address many challenges of modem astron- 
omy, but more detailed investigations (e.g. with specific 
chemical composition) based on synthetic spectra require 
a local implementation of synthetic spectroscopy codes, 
that are available off-the-shelf. The input data (line lists) 
is not yet consistently included in the VO, and the data 
model is still under development. The situation is better 
for atomic data, with the existence of the VALD database. 
Molecular data is scattered on the web and the best way 
so far to assemble line lists is to contact people that have 
done it for their own work. Some challenges are ahead, 
like the automatic classification of spectra, and the deter- 
mination of stellar parameters on samples of billions of 
spectra. New methods are being tested that may allow 
to do this reliably on a large scale. These tools could be 
included in the VO. Other tools that should be included 
in the VO are correlation and convolution (with rota- 
tion, macroturbulence, and instrumental profiles) tools 
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that could be used on both observed and calculated spec- 
tra. And finally, n all instances one has to make sure that 
proper credit is given to all developers of such tools, mod- 
els and calculations, if one wants to keep aUve a spirit of 
sharing resources in the community. Credit is unfortu- 
nately needed to get jobs and money. An possible solu- 
tion was given at the conference: forward the informa- 
tion on developers in the metadata sent to users, and in- 
form journal editors and referees that they should ensure 
proper referencing of the source of data extracted from 
the VO. And finally, it is not easy to go mining into cal- 
culated (or observed data), and we will have to keep on 
learning to use tremendous amounts of data in an edu- 
cated way. 
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